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Abstract 

New spectra from the FINUDA experiment of the Non-Mesonic Weak Decay (NMWD) proton 
kinetic energy for ;\Be, ^B, a 2 C, a. 3 C, A 5 N and A 6 are presented and discussed along with the 
published data on A He and A Li. 

Exploiting the large mass number range and the low energy threshold (15 MeV) for the 
proton detection of FINUDA, an evaluation of both Final State Interactions (FSI) and the 
two-nucleon induced NMWD contributions to the decay process has been done. Based on this 
evaluation, a linear dependence of FSI on the hypernuclear mass number A is found and for the 
two-nucleon stimulated decay rate the experimental value of r2/r p =0.43±0.25 is determined 
for the first time. A value for the two-nucleon stimulated decay rate to the total decay rate 
r2/rNMWD=0.24±0.10 is also extracted. 
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1. Introduction 



The Non-Mesonic Weak Decay (NMWD) of A-hypernuclei has been studied quite 
extensively both experimentally and theoretically since the early days of hypernuclear 
physics [1,2]. If, on one hand, studies on the structure of hypernuclei provide interesting 
information on the hyperon-nucleon strong interaction [3], on the other hand the weak 
decay of hypernuclei represents the easiest way to get information on the baryon-baryon 
strangeness-changing weak interactions [4,5]. 

In A-hypernuclei the NMWD occurs through processes which involve a weak inter- 
action of the A with one or more nucleons. Sticking to the hadronic vertex A — > irN, 
if the emitted pion is virtual it can be absorbed by the nuclear medium, resulting in a 
non-mesonic decay of the following types: 

An nn (T n ) , (1) 

Ap -> np (r p ) , (2) 

ANN ->■ nNN (T 2 ) . (3) 

The total weak decay rate of a A-hypernucleus is then: 

Tt = Tmwd + Tnmwd , (4) 

where rMWD=r T - + T^o is the total mesonic weak decay rate and rNMWD=ri + ^2, 
with Ti= r„ + r p ; the hypernucleus lifetime is r = H/Tt- 

The channel (3), referred to as two-nucleon induced decay and suggested in Ref. [6], 
can be interpreted by assuming that the pion from the weak vertex is absorbed by a pair 
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of nucleons, correlated by the strong interaction. Note that the non-mesonic processes 
can also be mediated by the exchange of mesons more massive than the pion. 

Several important experimental advances in NMWD study have been made in recent 
years, which have established more precise values of the neutron- and proton-induced 
decay rates T n and T p , solving in this way the long-standing puzzle on the T n /T p ratio 
[7]- ' 

Despite this progress, one should observe that no experimental evidence has been 
obtained yet of the two-nucleon stimulated decay, with the exception of two indirect 
results [8,9]. Various theoretical papers were dedicated to the calculation of the rates 
and the nucleon spectra for this mechanism [6,10-16], predicting a significant contri- 
bution to NMWD. The first calculation of the two-nucleon induced decay rates was 
performed within a nuclear matter scheme introducing a phenomenological description 
of the two-particle two-hole polarization propagator [6]. This model was improved in 
Rcfs. [10,11], where finite nuclei were considered in the local density approximation, and 
then a microscopic approach was developed [13-16]. In Refs. [6,10-12] only the decay 
channel Anp — > nnp was considered, while in the microscopic approach of Refs. [13-16] 
all three two-nucleon induced channels, including Ann — > nnn and App — > npp, were 
taken into account. In particular, a very recent contribution [16] studied the effects of 
Pauli exchange terms in the two-nucleon stimulated NMWD, thus giving more reliable 
evaluations of Fnmwd and I^/Tnmwd- 

The FINUDA experiment performed a complete analysis of the proton energy spectra 
following the NMWD of ^He, 7 A U, 9 A Be, ^B, \ 2 C, \ 3 C, \ 5 N and \ 6 hypernuclei. Results 
for ^He, ^Li and \ 2 C have already been published [17]. Exploiting the large systematics, 
we present here the results of a model-independent analysis of the contributions of Final 
State Interaction (FSI) and two-nucleon induced decays to the nucleon spectra. This 
study leads to the experimental indication of the two-nucleon induced decay process. 

2. The experimental and analysis method 

The results reported in the present Letter have been obtained by analyzing the data 
collected by the FINUDA experiment from 2003 to 2007 and correspond to an integrated 
luminosity of ~ 1.2 fb" 1 . The experimental method is briefly described here, while fur- 
ther details are reported in Ref. [17]. 

FINUDA is a fixed target experiment installed at one of the two interaction points of the 
DA$NE e+e~ ^-factory of Laboratori Nazionali di Frascati (INFN-Italy) . A detailed 
description of the experimental apparatus can be found in Rcfs. [18,19]. 
The FINUDA spectrometer has a cylindrical geometry and consists of several sensi- 
tive layers which are used for particle localization and identification (dE/dx). From 
the beam axis outwards, three main parts can be distinguished, as described below. 

- The interaction /target region: it is located at the center of the apparatus. Here the 
kaons coming from the 4> decay are identified by means of a barrel of 12 thin scintillator 
slabs (TOFINO), surrounded by an octagonal array of Si microstrips (ISIM) facing 
eight target modules. 

- The external tracking device: it consists of four layers of position sensitive detectors 
arranged in a coaxial geometry and immersed in a He atmosphere to minimize the 
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multiple scattering. A decagonal array of Si microstrips (OSIM) faces the target tiles; 

the next tracking devices are two octagonal arrays of low mass drift chambers (LMDC) 

and a stereo system of straw tubes (ST). 
- The external time-of- flight detector: it is a barrel of 72 scintillator slabs (TOFONE) 

providing signals for the trigger and time-of-fiight measurements. 
The whole apparatus is placed inside a uniform 1.0 T solenoidal magnetic field. 
The present analysis was performed on events collected out of three 12 C targets (1.7 mm 
thick, mean density 2.265 g cm -3 ), two 6 Li ones (90% enriched, 4 mm thick), two 7 Li 
ones (4 mm thick), two 9 Be ones (2 mm thick), a 13 C one (99% enriched powder, 10 mm 
thick, mean density 0.350 g cm -3 ) and a D2O filled one (mylar walled, 3 mm thick). 
The study of the proton spectra following the NMWD was performed by selecting all the 
events with two particles emitted in coincidence: a n~ carrying the information of the 
hypernuclcus formation (whether in its ground state or in a low lying excited one) and 
a proton coming from the same K~ interaction vertex which gives the signature of the 
NMWD. 

Negative pions and protons were identified by means of both the specific energy loss in 
OSIM and in the LMDC's and, if present, by the mass identification from the time-of- 
flight system (TOFINO-TOFONE). 

The 7r - 's were selected by requiring tracks identified by four hits (long tracks) , one in 
each of the FINUDA tracking detectors, whereas for proton identification we required 
not only long tracks but also tracks which do not reach the ST system (hereby indicated, 
for brevity, as short tracks). The detection threshold on the energy of protons was 15 
MeV, a lower limit never reached by previous experiments. 

By using this multiple particle identification and the described selection criteria, a 7r~ 
momentum resolution of Ap/p ~ 1% FWHM in the momentum region 260-280 MeV/c 
was achieved. The measured proton identification efficiency was 90% and the proton 
energy resolution was AE/E ~ 2% FWHM at 80 MeV. 

The data analysis starts with the selection of the ir~ momentum region corresponding to 
the A-hypernucleus formation in its ground state or in a low lying excited one decaying to 
the ground state by electromagnetic emission. The details of the intervals determination 
are exhaustively described in Ref. [20]. The momentum and binding energy intervals 
selected to identify the formation of the different hypernuclei under study are reported 
in Table 1. In order to increase the statistics for ^ 2 C the 7r~ momentum selection was 
extended to all the A bound region instead of including only the ground state [17]. 

The proton spectra were obtained by selecting all the protons detected in coincidence 
with a 7T _ in the momentum region reported in the third column of Table 1. 
In the raw ir~ and proton spectra background is mainly due to the K~(np) — > 
absorption on two correlated nucleons in the absorbing nucleus, followed by the in-flight 
£~ — > n~n decay. Out of the several processes in K~ absorption at rest that lead to a 
continuous spectrum of w~ , this is the only one affecting the region of the bound states of 
hypernuclei. This background was simulated, reconstructed with the same reconstruction 
program used for the real data and subtracted from the raw proton spectra. 
Another background source for ir~ and proton coincidence spectra is the A q f production 
and decay. However, the protons coming the in-flight A q f decay do not contribute to the 
proton spectra coming from the NMWD because they are below the FINUDA proton 
detection threshold. 

A crucial point of the analysis is then the acceptance correction of proton spectra; the 
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target 


hypernucleus 


(MeV/c) 


B.E. 

(MeV) 


references 


6 Li 




272 - 278 


0.63-5.99 


[21,22] 


7 Li 


5Hc 


267 - 273 


-3.02 - 1.85 




7 Li 


iLi 


273 - 279 


1.85 - 7.45 


[23-25] 


9 Bc 


«Be 


280 - 286 


1.50 - 7.00 


[23,24,26] 


12 C 




258 - 264 


-2.00 - 2.75 


[18,27] 


12 c 


A C 


264 - 273 


2.75 - 14.00 


[18,27] 


13 C 


i 3 c 


274 - 292 


-2.37 - 14.30 


[24] 


16 




265 - 270.5 


0.0 - 4.90 


[24] 


16 


ifo 


270.5 - 282 


4.90 - 15.40 


[24] 



Table 1 

Summary of the momentum and binding energy intervals selected to identify the formation of various 
hypcrnuclcar systems. First column: target nucleus; second column: formed hypernucleus; third column: 
momentum interval; fourth column: binding energy interval; fifth column: references of previous missing 
mass spectroscopy experiments. 

acceptance, e, was determined by taking into account the apparatus geometry, the effi- 
ciency of the FINUDA pattern recognition algorythm, the trigger request and the quality 
cuts applied to real data. 

The acceptance function, R = 1/e, for protons was evaluated target by target in different 
apparatus sector and for the energy it has a negative quadratic exponential behaviour 
in the 0-50 MeV range and flattens above 50 MeV. A complete account of the analysis 
procedure for ^He, ^Li and \ 2 C is given in Ref. [17]. 

3. Results and Discussion 

In Fig. 1 the kinetic energy spectra of protons coming from the NMWD of ^He, ^Li, 
^Be, J^B, ^ 2 C, ^ 3 C, \ 5 N and ^ 6 are shown. All the spectra are background subtracted 
and acceptance corrected; the errors in the spectra are statistical only and include both 
the contributions from background subtraction and acceptance correction. We estimated 
a systematic error of less than 5%. 

The analysis of data collected on several medium mass hypernuclei is a suitable way to 
extract more information on the NMWD process and is the relevant point of the present 
Letter. 

Even though affected by considerable errors, in particular in the low energy region, the 
proton kinetic energy spectra show a clear trend as a function of the hypernuclear mass 
number A (from 5 to 16): a peak around 80 MeV (which corresponds to about half of the 
Q-value for the free Ap — > np reaction) is broadened by the Fermi motion of nucleons 
and more and more blurred as A increases. The peak is smeared, on its low energy side, 
by a rise that can be ascribed to FSI and two-nucleon induced weak decays. 
Each proton spectrum from 80 MeV onwards was fitted to a Gaussian function. The 
mean values of the Gaussian are reported in Table 2 for each hypernucleus. The FWHM 
of the Gaussian fits is different for each spectrum and this variation is not related to our 
energy resolution for protons but it is due to the nucleon Fermi motion. 
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iZ^n + p+ A - 2 (Z-l) , 



(5) 



one has to specify the binding energy of the residual nuclear state A ~ 2 (Z — 1). A few 
different cases can occur: a single residual nucleus, in the ground state or in an excited 
state, or a more complicated system composed of two or more nuclear fragments (possibly 
in excited states) can be produced by the decay. To our knowledge there are no detailed 
calculations on this subject. We have verified that, in the two cases of ^He and ^Li, 
the mean values obtained by the fit are close to those obtained by assuming that the 
residual nucleus A ~ 2 (Z — 1) is emitted in the ground state ( 3 H and 5 He respectively). The 
calculated values in such cases are 70 MeV and 79 MeV, respectively. With increasing A, 
the fitted peak centers become increasingly smaller than the values calculated by means 
of a three body phase space calculation (more probable value) . This circumstance seems 
to indicate that decays, in which the final nuclear state consists of two or more fragments, 
are more probable for larger A. 



Hypcrnuclcus 


mean value 




68.5 ± 4.1 


iLi 


76.7±5.2 


\Be 


78.2±6.2 


A X B 


75.1±5.0 


12c 
A ° 


80.2±2.1 


i 3 c 


83.9±12.8 


A 5N 


88.1±6.2 


A 6 © 


93.1±6.2 



Table 2 

Mean values of the Gaussian fit performed on the proton energy spectra following the NMWD of ^He, 
iLi, \Be, iiB, fC, 13C, A 5 N and 

Let us now consider the low energy rise in the proton spectra, which we consider to be 
due to the combined contributions of FSI and two-nuclcon induced NMWD. If it were due 
to FSI effects, one should naturally expect that the structure at 80 MeV should still exist, 
but reduced as number of events since a fraction of them would be degraded in energy 
and contribute to the low energy part of the spectra. By straightforward arguments we 
may assume that the number of protons suffering a FSI should increase linearly with A 
and then the related amount of events in the tail due to FSI. 

The additional mechanism that could explain the low energy behaviour is a contribution 
of the two-nucleon induced NMWD; if the energy released in the weak decay is shared 
by three nucleons, a low energy rise may show up even for hypernuclei as light as %He. 
This contribution to the NMWD should be assumed, as first hypothesis, as independent 
on (or slowly dependent on) A. Theoretical calculations [4] support such a constancy of 
r2 /Tnmwd for a wide hypernuclear mass range. 

A way of disentangling between the two mechanisms is that of studying the behaviour 
of the low energy tail as a fuction of A as we will describe in the following. 

To estimate the contribution of two-nucleon induced NMWD we first assume that 
the part of proton spectrum beyond the peak mean value at about 80 MeV is due to 
protons coming from the Ap — > np reaction and that the two-nucleon induced channel 
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can be neglected in this region. This assumption is confirmed by theoretical calculations, 
which estimated that beyond 70 MeV the protons from the two-nucleon stimulated decay 
are about 5% of the total number of two-nucleon induced decay protons [12]. On the 
contrary, the part of spectrum below the peak mean value is due to protons from both 
Ap — > np and two-nucleon induced decays. 

As pointed out before FSI also affects the spectra. Let us then denote with ]>}F SI - low 
(N FSI ~ hlgh ) the difference between the number of protons detected in the region below 
(beyond) the peak of the spectrum and the (unknown) number of primary protons which 
originates from one- and two-nucleon induced weak decays in the same region. Available 
calculations [12] predict that N FSI_low is a positive quantity while N FSI ~ hlgh is negative: 
FSI tends to remove protons from the high energy part of the spectrum while filling the 
low energy region. 

According to the above assumptions and definitions, for the total number of protons 
Aj ow in the low energy region (below the center of the peak) and Ahi g h in the high energy 
region (beyond the center of the peak; blue filled areas in Fig. 1) of the spectra we thus 
have: 

A low = N(Ap -> np)/2 + N(Anp -> nnp) + A^ FSI - low , (6) 
A high = N(Ap -> np)/2 + ATfSi-high _ (7) 

where N(Ap — > np) is the total number of Ap — > np decays (which we assume to be evenly 
distributed in the Gaussian peak, i.e., half in Ai ow and half in Ahigh) and N(Anp — > nnp) 
is the total number of Anp — > nnp decays (which occur only in Ai ow ) contributing to the 
spectra. One can also write: 



N(Anp — >■ nnp) T np r 2 
N(Ap np) ~ TV ~ IV 



(8) 



Note that in Eq. (6) and in the last equality of Eq . (8) the two-nucleon induced NMWD 
has been assumed to be dominated by the Anp — >■ nnp channel: r 2 = T(Anp — > nnp) + 
T(App — > npp) + T(Ann — > nnn) = T np + T pp + T nn ~ T np (the recent microscopical 
calculation of Rcf. [16] obtained T np : T pp : T nn = 0.83 : 0.12 : 0.04). 
Consider now the ratio: 

A low 0.57V(Ap ->■ np) + N(Anp ->■ nnp) + N^ sl - low 

R ~ Alow + Ahigh ~ N(Ap ^ np) + N(Anp ->■ nnp) + iV p FSI - low + iV p FSI - high ' ^ 

In Fig. 2 the experimental values of this ratio are plotted as a function of A. Due to 
FSI effects, both Ai ow and Ahigh given by Eqs. (6) and (7) arc expected to be propor- 
tional to A: one can thus perform a fit of the data of Fig. 2 with a function R(A) = 
(a' + b'A)/(c' + d'A). From this fit it turns out that d'/c' < 10~ 3 , thus a linear fit, 
R(A) = a + bA, can be considered. The result of such a fit (a = a'/c' = 0.654 ± 0.138, 
b = b'/c' = 0.009 ± 0.013 and x7Mf = 0.298/6) is shown in Fig. 2. 
Now, we note that we can write r 2 /r p = (r 2 /Ti)(l + r„/r p ). Since theory and experi- 
ment support the indipendence of r 2 /ri [4] and r„/r p [4,8] on A, we may assume that 
r2/F p is constant. 
So one can then rewrite Eq. (9) as: 
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0-5 + ^ 
R(A) = =±£ + bA 



(10) 




Fig. 2. The ratio 

^low/(^low ~t~ ^high) as a function of the hypernuclear mass number. 

Eq. (10) can be solved for T 2 /T p for each hypernucleus, getting: 

T 2 _ [R(A) - bA] - 0.5 
f£ ~ 1 - [R(A) - bA] ■ 



(11) 



The values obtained for r 2 /r p for our hypernuclei are all compatible with each other 
within errors. The final result can thus be given by the weighted average: 



0.43 ±0.25. (12) 



To determine r 2 /rNMWD we need to know the T n /T p ratio. Indeed: 

r 2 = r 2 /r p 
Tnmwd (r„/r p ) + 1 + (r 2 /r p ) 



(13) 



Theoretical calculations [7] predict r n /T p values between 0.3 and 0.5 for |He and \ 2 C 
and recent experimental results [8] give: T n /T p = (0.45 ± 0.11) for ^He and T n /T p = 
(0.51 ±0.13) for \ 2 C. To determine a unique value of T n /T p that could be reasonable for 
all the studied hypernuclei, we take the weighted average of the two cited experimental 
results, r„/r p = (0.48 ±0.08). Using this value together with our determination of r 2 /r p 
we obtain: 

0.24 ±0.10. (14) 



Tnmwd 



One may wonder that the above value (14) is affected by the experimental threshold 
of 15 MeV. However it was shown [14] that proton spectra from two-nucleon induced 
NMWD, having a three body like phase space distribution, are almost vanishing below 
15 MeV. On the contrary the proton spectra from FSI have a monotonic increase going 



9 



down in energy. Than we expect that the threshold at 15 MeV does not substantially 
affect the final value of I^/Tnmwd- The value derived for the I^/Tnmwd ratio supports 
both theoretical predictions [10-16] and the latest experimental results of Ref. [8,9], which 
claim a large contribution of the two-nucleon induced channel to the total NMWD rate. 
However, the value obtained for I^/Fnmwd suggests a smaller contribution than what 
quoted by Bhang et al. [8], that for ^ 2 C reported a contribution of the two-nucleon 
induced process as large as 40% of the total NMWD width. We remark that this value 
was inferred from the fit of a KEK proton spectrum measured with good statistics but 
poor energy resolution. The KEK spectrum is quite different from the one reported here. 
Experimental reasons for such a discrepancy are discussed in Ref. [17]. 

4. Conclusions 

We reported a systematic study, performed by the FINUDA experiment, of the proton 
kinetic energy spectra following the NMWD of A-hypernuclei in the A = 5 -f- 16 mass 
range. 

Taking advantage of the low kinetic energy threshold for proton detection (15 MeV) and 
of the large mass number range of the studied hypernuclei, an evaluation of the nucleon 
FSI effects and of the two-nucleon induced decay contribution to the decay process was 
performed. 

A linear dependence of FSI on the hypernuclear mass number A was found and an 
experimental value of the ratio r 2 /r p = (0.43 ± 0.25) was determined for the first time. 
This result provides a clear experimental indication of the relevant role played by the 
two-nucleon induced mode in the NMWD of hypernuclei: this channel contributes to 
about 24% of the Tnmwd for the studied hypernuclei. 
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